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We report low-temperature carrier transport properties of a series of nominally uncompensated 
neutron-transmutation-doped 70 Ge:Ga samples very close to the critical concentration JV C for the 
metal-insulator transition. The nine samples closest to N c have Ga concentrations TV in the range 
0.99-/V c < N < 1.01iV c . The electrical conductivity a has been measured in the temperature range 
T = 0.02 - 1 K. On the metallic side of the transition the standard a(T) = a + bT q with q = 1/2 
was observed for all the samples except for the two that are closest to N c with N between N c and 
1.00157V C . These samples clearly show q — 1/3. 

An extrapolation technique has been developed in order to obtain the zero-temperature conduc- 
tivity a(0) from o(T) with different dependence on T. Based on the analysis, v ~ 0.5 in the familiar 
form of cr(0) oc [N/N c — 1)" has been found. On the insulating side of the transition, variable 
range hopping resistivity p(T) oc exp(Tb/T) p with p = 1/2 has been observed for all the samples 
having N < 0.991iV c . In this regime T oc (1 - N/N c ) a with a w 1 as N -> N c . The values of 
To agree very well with theoretical estimates based on the modified Efros and Shklovskii relation 
fcjjTo ~ (2.8e 2 /47re Ko£o)(l - N/N c ) a , where k and £ 

are the dielectric constant and the Bohr 
radius, respectively. The insulating samples very close to the transition (0.9917V C < N < N c ) exhibit 
quite a different behavior. In this range 1/p increases rapidly as TV changes from 0.991iV c to N c . 
The relevance of our findings to the collapsing of the Coulomb gap is discussed. 

PACS numbers: 71.30.+h, 72.80.Cw 



I. INTRODUCTION 

The daping-induced metal-insulator (MI) transition in semiconductors has been studied extensively in the past few 
decades. Ertil However, there still remain a number of major theoretical and experimental challenges. Measurements 
of the electrical conductivity cr(T) as a function of temperature near the MI transition are fundamental to the 
understanding of the roles of potential disorder and electron-electron interaction. The zero-temperature conductivity 
tr(0) obtained from an appropriate extrapolation of the temperature dependent cr(T) to zero temperature is evaluated 
as a function of doping concentration N immediately above the critical concentration N c for the MI transition; 

a(0)=a (N/N c -iy, (1) 

where ctq is the prefactor and v is the critical exponent. In_several strongly disordered systems, e.g., compensated 
single crystalline semiconductors [Ge:SbJj Si:P,B,ElGe:Ga,As,El Alo.3Gao.7As (Ref. |[)] and amorphous alloysjjtj v « 1 
has been found. These results are in good agreement with the predictionta for the transition driven by disorder. It was 
also found that in compensated Alo.3GaorfAs the dielectric constant on the insulating side diverges with the critical 
exponent of s ~ 2.3 near the transition,El i.e., s w 2v predicteo£3 for the disorder-induced transition holds. Thus 
there is strong evidence that the effect of disorder rather than electron-electron interaction plays the key role in the 
MI transition of compensated semiconductors. On the other hand, a critical, exp onent of-t' w 0.5-has been obtained 
with a number of nominally uncompensated semiconductors [Si:P,l 13 H 14 l Si:Aa Ma Ge:As,tZI Si:B,t3 Ge:Ga (Ref. ||)]. 
This value of v w 0.5 is significantly smaller than v w 1 — 1.3-4)redictedaEJ by the transition purely driven by 
the disorder. It also does not satisfy Chayes et aVs inequality^ v > 2/3 for transitions due to both disorder and 
electron-electron, interaction. In response to these discrepancies, several theoretical ideas supporting v w 0.5 have 
been proposed EiJ~E3 However, general agreement between the experimental results and theory has yet to be achieved 
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by any of the models H3"Ea The interesting observation reported commonly on uncompensated systems is the relatively 
wide range of N above N c (typically up to ln^Vc or larger) in which <r(0) can be fitted very well with Eq. (Q) with 
v w 0.5. Based on this observation, Fritzsche£3 proposed ajsiodel composed of one main transition accompanied by 
two satellite transitions, one on each side of N c . Stupp et allEi questioned the large critical region and found a narrow 
regime N c < N < l.liV c in which <j(0) of uncompensated Si:P exhibits v « 1.3. More recently, v » 1 was claimed also 
for uncompensated Ge:As by Shlimak et alt3 This recent trend of v moving from « 0.5 towl- 1.3 was ended by op; 
work on homogeneously doped, nominally uncompensated Ge:Ga, in which v ss 0.5 was established unambiguously.1121 
The exponents v w 1 — 1.3 claimed for melt-doped Si:P (Ref. [3(]) and Ge:As (Ref. |3l]) should be interpreted with 
great caution for the reasons we give in the following paragraphs. 

In the experiment reported here we probe the low-temperature electrical properties of nominally uncompensated 
Ge:Ga in the region extremely close to the MJ-,transition; 0.99^ < N < 1.017V C . This concentration regime has 
not been fully investigated in our earlier work.E^LFpi:. the . case of melt- (or metallurgically) doped samples that have 
been employed in most of the previous studies ,Ej~t§E3Ell the spatial fluctuation of N due to dopant striations and 
segregation can easily_be on the order of 1% across a typical sample for the four-point resistance measurement (length 
of ^5 mm or larger). l3 For this reason it is not meaningful to discuss physical properties in this truly critical region 
(e.g., \N/N C — 1| < 1%) based on the data taken with melt-doped samples. 

A precise determination of N in a melt-doped sample is also difficult due to the spatial fluctuation of N as well as to 
the limited accuracy of the existing method to measure N near the transition. The determination of N by Hall effect 
may be inaccurate due to the possible divergence of the Hall coefficient from unity near the transition. Resistivity 
measurements at two temperatures (4.2 K and 300 K) (Ref. [if]) to find N require an accurate calibration that cannot 
be established easily. I— . 

All Ge:Ga samples used in this work (and in our earlier studytLij) were prepared by neutron-transmutation doping 
(NTD) of isotopically enriched 70 Ge single-cwstals. Our NTD method inherently guarantees the random distribution 
of the dopants down to the atomic level.E3t£l The N for each sample is given by the thermal neutron fluence and 
its relation to TV has been accurately establishedta for 70 Ge. We prepared 13 new NTD 70 Ge:Ga samples with rune 
of them in the 0.99iV c < N < 1.01iV c region. The sample with the N closest to N c has N = 1.0004iV c . To our 
knowledge, neither experimental nor numerical studies on the MI transition have ever approached N c as close as this 
work has. Our study focuses on the analysis of the temperature dependence of <r(T) below 1 K on both sides of the 
transition; the insulating phase (N < N c ) and the metallic phase (N > N c ). We investigate the universality of the 
<t(T) in the metallic phase by introducing a numerical procedure. A quantitative discussion of cr(T) in the insulating 
phase will be given in the context of the variable range hopping conduction model. 

II. EXPERIMENT 
A. Sample preparation and characterization 

We first describe the preparation of the neutron-transmutation-doped 70 Ge:Ga samples for the low-temperature 
conductivity measurements in the critical regime of the MI transition. We use NTD finee it is known to produce 
the most homogeneous, perfectly random dopant distribution down to the atomic level.&tiS The Czochralski grown, 
chemically very pure 70 Ge crystal has isotopic composition [ 70 Ge]=96.2 at. % and [ 72 Ge]=3.8 at. %. The as-grown 
crystal is free of dislocations, p type with an electrically active net-impurity concentration less than 5xl0 n cm -3 . 
The thermal neutron irradiation leading to NTD was performed at the University of Missouri Research Reactor with 
the thermal to fast neutron ratio of ~30:1. Upon capturing a thermal neutron 70 Ge becomes 71 Ge which decays with 
a half-life of 11.2 days via electron capture to a 71 Ga acceptor. The small fraction of 72 Ge becomes 73 Ge which is 
stable, i.e., no further acceptors or donors are introduced. The post NTD rapid-thermal annealing at 650 °C for 10 sec 
removed most of the irradiation-induced defects from the samples. The short annealing time is important in order 
to avoid the redistribution and/or clustering of the uniformly dispersed 71 Ga acceptors. The concentration of the 
electrically active radiation defects-measured with deep level transient spectrometry (DLTS) after .the annealing is less 
than 0.1% of the Ga concetrationj 3 ^ i.e., the compensation ratio of the samples is less than O.OOI.lD The dimension of 
most samples used for conductivity measurements was 6x0.9x0.7 mm 3 . Four strips of boron-ion-implanted regions 
on a 6x0.9 mm 2 face of each sample were coated with 200 nm Pd and 400 nm Au pads using a sputtering technique. 
Annealing at 300 °C for one hour activated the implanted boron and removed the stress in the metal films. 

The Ga concentration N in our 70 Ge samples after NTD is given precisely by 

[ 71 Ga] (cm" 3 ) = 0.1155 x n (cm" 2 ), (2) 
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where n is the thermal neutron fluence.EZl The main goal of this study was to fill the gap in N between 1.840 and 
1.861xl0 17 cm -3 that was missing in our earlier workjia i.e., a precise control of n of the order of 0.1% is needed. 



Although obtaining such a precision in n seems dif 
pared the insulating samples for our previous study. 



,ult, we successfully used the following approach. When we pre- 
I we doped three 2 cm diameter wafers to N = 1.733xl0 17 cm~ 3 . 
For the present study 13 pieces were cut from two of the N = 1.733xl0 17 cm~ 3 wafers. Each of the 13 pieces were 
then irradiated a second time to cover the range N = 1.840 — 1.861 xlO 17 cm -3 with a neutron fluence resolution of 
n = 2.2xl0 15 cnr 2 which corresponds to N — 2.5xl0 14 cm -3 according to Eq. (||). 



B. Measurements 



The electrical conductivity measurements were carried out down to temperatures of 20 mK using a 3 He- 4 He dilution 
refrigerator. All the electrical leads were low-pass filtered at the top of the cryostat. The sample was fixed in the mixing 
chamber and a ruthenium oxide thermometer [Scientific Instrument (SI), RO600A, 1.4x1.3x0.5 mm 3 ] was placed 
close to the sample. To measure the resistance of the thermometer, we used an ac resistance bridge (RV-Elekroniikka, 
AVS-47). The thermometer was calibrated against 2Ce(N0 3 ) 3 • 3Mg(N0 3 ) 2 -24H 2 (CMN) susceptibility and against 
the resistance of a canned ruthenium oxide thermometer (SI, RO600A2) which was calibrated commercially over a 
temperature range from 50 mK to 20 K. We employed an ac method at 21.0 Hz to measure the resistance of the 
sample. The power dissipation was kept below 10~ 14 W, which is small enough to avoid overheating of the samples. 
The output voltage of the sample was detected by a lock- in amplifier (EG&G Princeton Applied Research, 124A). 
All the analog instruments as well as the cryostat were placed inside a shielded room. The output of the instruments 
was detected by digital voltmeters placed outside the shielded room. All the electrical leads into the shielded room 
were low-pass filtered. The output of the voltmeters was read by a personal computer via GP-IB interface connected 
through an optical fiber. 



III. RESULTS AND DISCUSSION 



A. Electrical transport and the critical conductivity exponent in the metallic samples 

The temperature dependence of the electrical conductivity mostly for the metallic samples is shown in Fig. |l[ The 
solid symbols denote the data taken with the samples prepared in this work and the open ones are the data reevaluated 
with most of the samples described in Ref. [lj]. Note that several samples are doped successfully in the immediate 
vicinity of N c . Mott's minimum metallic conductivity er m i n for Ge:Ga is estimated to be 7 S/cm using the relation 

CT m in = CM(e 2 /h)Ny a with Cm ~ 1/20 as it was used for Si:P.EB0 Figure [l] clearly shows that o of some of the newly 
prepared metallic samples takes values less than <r m i n even at finite temperatures. The critical exponent v in Eq. (|l|) 
is defined for the critical region [N/N c -1) < 1 through the conductivity at zero temperature cr(0). Experimentally, 
however, it is impossible to reach T = and a suitable extrapolation is required. 

The temperature variation of the conductivity is governed mainly by the electron-electron interaction and can be 
written as 

Act(T) = cr(T) - cr(0) = mVT, (3) 

where 

m = A/VB . (4) 

Here, A is a temperature independent constant and D is the diffusion constant, which is related to the conductivity 
via the Einstein relation 

a = (dn/dfi)e 2 D, (5) 

where (dn/dfj.) is the density of states at the Fermi level. In the limit of Aer(T) <C <r(0) ~ &(T), D can be considered as 
a constant, i.e., m is constant. Usually <r(0) is obtained by extrapolating <r(T) to T = assuming VT dependence based 
on Eq. (||). Such an analysis was performed in our earlier work since Acx(T) cx \fT was found for all the samples.E2l 
It should be pointed out, however, that the above inequality Aer(T) <C <r(0) is no longer valid as N approaches N c 
from the metallic side since <r(0) also approaches zero. In such cases m in Eq. (||) is not temperature independent and 
Aa(T) may exhibit a temperature dependence different from y/T. To examine this point for our experimental results, 
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we go back to Fig. |l| It is seen here that the Act(T) of the bottom five curves are not proportional to VT while 
Aer(T) of the other higher TV samples are well described by oc VT. The close-ups of cr(T) for the six samples with 
positive da/dT in the scale of VT and T 1 / 3 are shown in Figs. ||(a) and||(b), respectively. The upper and lower dotted 
lines represent the best fit using the data between 0.05 K and 0.5 K for the samples with TV = 1.912 x 10 17 cm~ 3 
and N — 1.861 x 10 17 cm~ 3 , respectively. Each fit is shifted downward slightly for easier comparison. From this 
comparison it is clear that a T 1 / 3 dependence holds for samples in the very vicinity of the MI transition. The opposite 
is true for the curve at the top. This means that the VT dependence in Eq. (^) is replaced by a T 1 / 3 dependence as 
the MI transition is approached. ._. 

A T 1 / 3 dependence close to the MI transition was predicted originally by Al'tshuler and AronovH They considered 
an interacting electron system with paramagnetic impurities, for which they obtained a single parameter scaling 
equation. At finite temperatures, they assumed a scaling form for conductivity according to the scaling hypothesis; 

o = ^ m/Lr), (6) 



where £ is the correlation length and Lx = V^D/ksT is the thermal diffusion length. When ^> £, /(£/£t) = 
A + B(£/Lt), which is equivalent to Eq. (0). In the critical region, where Lt <C £ — > oo, Eq. (0) should be reduced to 
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<* = C-j- ■ (7) 

Combining this equation and Eq. (^|)_£hey obtained a oc T 1 / 3 . More recently, the T 1 / 3 dependence has been predicted 
based solely on the effect of disordered and on the quantum interference £3 

Although the origin for the T 1 / 3 dependence in the present system near N c is unknown at this point, it is impor- 
tant that we find a method that allows the determination of cr(0) even when the temperature dependence of 
changes from VT to T 1 / 3 as N approaches N c . For this purpose we follow Al'tshuler and Aronov's manipulation!^ of 
eliminating m and D in Eqs. (ff)-® and obtain 



a(T) = o-(O) + m'VT/a(T) , (8) 



where ml — Ae^J (dn/dfi) , which is temperature independent. In the limit of Aer(T) <C <r(0) ~ c(T') , this equation 
gives the same value of cr(0) as Eq. (§) does. When cr(0) < a(T), it yields a T 1 / 3 dependence for a(T). Thus it 
is applicable to both VT and T 1 / 3 dependent conductivity. From today's theoretical understanding of the problem, 
Eqs. (|) and (§) are valid only fovr-Lx ^ £, and their applicability to the critical region is not ri*ar, because the 
higher-order terms of the f3 functionta which were once erroneously believed to be zero do not vanish.ESca Nevertheless, 
we expect Eq. to be a good expression for describing the temperature dependence of all metallic samples because it 
expresses both the VT and the T 1 / 3 dependences as limiting forms. Then, based on Eq. (||) we plot a(T) vs y/TJa(T) 
for the four close to N c samples in Fig. ||. As we see, the data points align on straight lines very well, which supports 
the adequacy of Eq. (|§|). The zero-temperature conductivity cr(Q) is obtained by extrapolating to T = 0. The curve on 
the top of Fig. ||is for the sample with the lowest N among the ones showing VT dependence at low temperatures, i.e., 
this sample has the largest value of A<r(T)/cr(0) among VT samples. The value of cr(0) obtained for this particular 
sample using Eq. (^) differs only by 0.6% from the value determined by the conventional extrapolation assuming 
Eq. (|^). This small difference is comparable to the error arising from the choice of the temperature range in which 
the fitting is performed. Therefore the new extrapolation method proposed here is compatible with the conventional 
method based on-the oc VT extrapolation. A different method for the determination of er(0) was proposed recently 
by Shlimak et alxA but it requires u{T) — a + bT 1 / 3 for all samples with exactly the same b. Such a strict condition 
is not met in Ge:Ga and in many other systems as we will show later in Fig. |5|. 

Based on our new analysis, the MI transition was found to occur between the first and second samples from the 
bottom in Fig. || Here N c is fixed already within an accuracy of 0.16% corresponding. .to -the fractional difference in 
N between the first and second samples from the bottom, i.e., unlike the case for Si:P,E3c3'Cj the determination of the 
critical conductivity exponent will not be affected by the ambiguity in the value of N c . Figure [| shows the <r(0) as a 
function of N/N c - 1 with an excellent fit by Eq. @ (dotted line) with v = 0.50 ± 0.04 and N c = 1.860 x 10 17 cm^ 3 
all the way down to (N/N c — !,)-■= 4 x 10~ 4 . The clear demonstration of the same v rj 0.50 in our previous workEa 
was criticized by Shlimak et alsEi for our doping level's not being close enough to N c . The present work shows that 
the critical exponent is indeed ~ 0.5 for nominally uncompensated Ge:Ga. We note that v — 0.46 ± 0.18 w 0.5 is 
obtained even when fitting only the results obtained with the four samples closest to the transition. As was mentioned 
in the Introduction, Eq. (|l|) with v w 0.5 holds for many nominally uncompensated crystalline semiconductors for 
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a relatively wide range of N above N c . The cr(0) of Ge:Ga shown in Fig. g can be fitted very well with a single 
exponent v w 0.50 over three orders of magnitude in N/N c — 1. In fact it was shown in Ref. |l9| that ^ w 0.5 holds up 
to N = IAN C . 

In order to compare the low-temperature transport properties of Ge:Ga with other systems, we evaluate the con- 
centration N* where the sign of da/dT changes. N* of our p-type Ge:Ga lies between 1.04A^ C and 1.08./V c _syhile 
larger values of 1.2iV c < N* < 1.3N C have been reported for n-type germanium; Ge:Sb (Ref. g|) and Ge:AsB The 
magnitude of Aa(T) in Ge:Ga is considerably smaller than that of Ge:Sb (Ref. ^) for samples with approximately 
the same N/N c — 1. A number of properties related to the band structure, e.g., the valley degeneracy, strength of 
the spin-orbit interaction, the degree of the intervalley scattering, etc., can change the low-temperature transport 
properties of doped semiconductors. The difference in the behavior of cr(T) at finite temperature between p- and 
n-type Ge may be understood in such contexts. Concerning the critical behavior of <r(0) at the MI transition, it is 
usually thought to depend on the universality class to which the system belongs, and can vary depending on the 
strength of the spin-orbit scattering or of the spin scattering. From the experiments so far done, including the present 
one on doped semiconductors except n-type Ge, we conclude, however, that the critical exponent v w 0.5 applies, 
irrespective of the systems as long as the compensation is not4*nportant. Regarding n-type Ge, v sa 1 was reported 
in Ge:As (Ref. [Jl]) and Ge:Sb (Ref. |3l| ) and v w 0.9 in Ge:Sb.c£l In order to verify whether this is truly the case, an 
investigation of n-type NTD 74 Ge:As is important. It is also interesting to point out that iV* of Ge:Ga is very similar 
to N* found in both p- and n-type Si. The N* for Si:B (Ref. g?]) is about 1.08JV C and for Si:P (Refs. |c] and g|) lies 
between 1.03iV c and 1.2N C . The coefficient m in Eq. (||) is compared in Fig. || for Ge:Ga, Si:P, and Si:B systems. 



B. Variable range hopping conduction in insulating samples 

The temperature dependence of the resistivity of insulating samples is shown in Fig. ^. Shklovskii and Efros have 
shown for insulating samples that a parabolic shaped energy gap (kuawn as the Coulomb gap) exists in the sigle- 
particle density of states in the immediate vicinitv-pf the Fermi level.c3 The variable range hopping resistivity p for 
the excitation within the Coulomb gap is given byc3 



where po is a prefactor, p = 1/2, and 



p = p exp[(T /r)P], (9) 



k B T^-^-^—. (10) 



47re k(N)£(N) 



k(N) and £(N) are the dielectric constant and localization length, respectively. Moreover, k(N) = kq(1 — N/N c ) s 
and £(N) = £o(l — N/N c )~^ as N approaches N c from the insulating side so that To becomes 

k B T a * 2 - 8e \ (1-N/N c r, (11) 
47re Koto 

where a = s + £ is to be determined experimentally in this study. Because the width of the Coulomb gap Acg 
depends also on N via k(N) as Acg oc [k(N)]~ 3 ^ 2 , it collapses rapidly as N approaches to from the insulating 
side. When Aqq becomes sufficiently small near N c , the excitation energy for hopping given by the thermal energy 
can become larger than Acg In this case the density of states may be considered to be constant around the Fermi 
level and the Mott variable range hopping with p = 1/4 in Eqt-|(|9|) is expected to be observed. Such a crossover from 
p=l/2top=l/4asiV approaches N c was observed in Si:P£il It is of great interest to see if such a crossover exists 
in our homogeneously doped Ge:Ga system. Figure shows the values of p found from the calculation of dine/ dlnT 
where e = — dlnp/ dlnT and/or from the direct fitting of curves shown in Fig. |^ by Eq. (||). p w 1/2 is obtained for 
the samples having N < 0.991^, i.e^they will be analyzed in the framework of Shklovskii and Efros's theory for 
the hopping within the Coulomb gap.CJ In Figs. |^(a) and |^(b), To and po, respectively, are platted as a function of 
(1 — N/N c ) for the samples with N < 0.991^. As was already shown in our previous workjij the best fit of T to 
Eq. ( |Tl| ) is obtained here with a « 1. Based on this finding, we calculate To using Eq. ( pd| ) with a = 1, kq = e = 16, 
and £o—f 47reoefi. 2 /m*e 2 = 8 nm, where e is the dielectric constant of Ge and m* is the effective mass of the electron 
in GeH The calculated (not fitted) T , which is shown by the dotted line in Fig. ||(a), agrees very well with the 
experimentally determined To, supporting the quantitative validity of the theoretical expression for Tq. 



In some of the earlier studies, the constant 2.8 in Ea^JlO) had to be adjusted to much smaller values in order 
to obtain an agreement with experimentally found Tp ,E§Hln Fig. H(b), po is shown as a function of 1 — N/N c . 
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The prefactor po shows no critical behavior and it approaches near N c a value very close to the inverse of Mott's 
minimum metallic conductivity denoted by the dotted line. Finally we turn our attention to p of the samples having 
0.9917V C < N < N c in Fig. 0. In this regime lying very close to N c , 1/p increases rapidly as N approaches N c due 
to the collapsing of the Coulomb gap. However, 1/p does not approach a constant value of 4 expected for the Mott 
variable range hopping conduction. In our analysis the temperature dependence of the prefactor po, which can be 
significant near N c , is neglected. Therefore further analysis taking into account the appropriate dependencies of po 
on T is important. Unfortunately we cannot perform such an analysis with— the accuracy needed at this point since 
the theoretical models proposed so far on po do not agree with one another O 



IV. CONCLUSION 



We have measured the electrical conductivity of nominally uncompensated neutron-transmutation-doped isotopi- 
cally enriched 70 Ge:Ga samples. Approaching the transition from the metallic side, we find that the temperature 
dependence of the form a(T) = a + bT q with q = 1/2 is replaced by q = 1/3. We introduce a method for finding 
<r(0) which is consistent with the conventional y/T extrapolation. The critical conductivity exponent v w 0.5 for 
p-type germanium has been fully confirmed. On the insulating side of the MI transition, the standard relation for 
the variable range hopping resistivity p(T) oc exp(To/T) p with p = 1/2 is observed for N < 0.991iV c . Shklovskii and 
Efros's expression for Tq agrees quantitatively with our experimentally found Tq. 
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FIG. 1. Electrical conductivity as a function of T 1 / 2 for NTD 70 Ge:Ga. From bottom to top in units of 10 17 cm -3 , 
the concentrations for the samples denoted by solid symbols are 1.853, 1.856, 1.858, 1.861, 1.863, 1.912, 2.210, and 2.232, 
respectively. Open symbols are the data taken on the samples used in our previous work (Ref. FTeJ) . 



FIG. 2. Conductivity as a function of (a) T 1 ^ 2 and (b) T 1 / 3 , respectively, near the MI transition. From bottom to top in 
units of 10 17 cm -3 , the concentrations are 1.853, 1.856, 1.858, 1.861, 1.863, and 1.912, respectively. The upper and lower dotted 
lines in each figure represent the best fit using the data between 0.05 K and 0.5 K for the first and the third curves from the 
top, respectively. Each fit is shifted downward slightly for easier comparison. 



FIG. 3. Conductivity a as a function of (T/<t) 1//2 . From bottom to top in units of 10 17 cm 3 , the concentrations are 1.858, 
1.861, 1.863, and 1.912, respectively. The solid lines denote the extrapolation for finding cr(0). 
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FIG. 4. Zero-temperature conductivity a(0) vs the dimensionless distance N/N c — ', 
logarithmic scale. The dotted line represents the best power-law fit by cr(0) oc (N/N c — 
symbols are from our previous work (Ref. O). 



1 from the critical point on a double 
1)" where v = 0.50 ± 0.04. The open 



FIG. 5. Coefficient m defined in Eq. (||) as a function of N/N c — 1; Ge:Ga of this work (•), Ge:Ga of the previous work 
(Ref. |l|) (o), Si:B (Ref. ^ (A), and Si:P (Ref. (□). 

FIG. 6. The logarithm of the resistivity as a function of T~ 1,/2 for insulating samples. The triangles denotes the data from 
Ref. [Is] with the fit by Eq. ^ (solid line). The samples of the present study are represented by circles and the concentrations 
from top to bottom in units of 10 17 cm" 3 are 1.840, 1.842, 1.843, 1.848, 1.850, 1.853, 1.856, and 1.858, respectively. 



FIG. 7. The inverse of the exponent p defined by Eq. fflh vs concentration. The open circles are from Ref. p_' 



FIG. 8. On a double logarithmic scale, To and po are plotted as functions of 1 — N/N c in (a) and (b), respectively. The open 
symbols are after Ref. |l9[ The dotted line in (a) is a calculated To using Eq. ( |l"l| ) with a = 1. The dotted line in (b) represents 
the inverse of the Mott minimum metallic conductivity. 
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